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ABSTRACT 


This paper presents a comprehensive overview of the analytical structure 
of sliding polynomials. Emphasis is placed on explanation so that researchers 
with limited background in numerical analysis can apply the methods. 
Presentation is generalized to assist readers from various disciplines in 
preparing data for interpretation. Polynomials of increasing complexity were 
previously presented in separate publications. The intent here is to develop 
an appreciation of the common elements underlying three levels of complexi- 
ty, starting with simple linear forms. Continuity of the first derivative of 


interpolation arcs is added first, folk 
second derivative. The three forms 
operators using least squares. 

INTRODUCTION 

Interpolation and smoothing are routinely re- 
quired to reduce and process experimental data. 
Choice of a particular method depends on the 
amount of data and on awareness of available 
methods. DuChateau et al. (1972) and Kimball 
(1974) dealt with interpolation and smoothing 
techniques that require substantial amounts of 
experimental data. Such data may be generated 
either in an extensive cyclic time domain or by 
repetition in a limited time or space domain. 
Savilzky and Golay (1964) discussed smoothing 
and differentiation of data using convoluting 
functions in weighted moving average tech- 
niques. Again, availability of substantial a- 
mounts of data is implied, generated by basically 
continuous physical experimentation. The author 
(1961, 1962, 1967a) showed how sliding polyno- 
mials could be used for interpolation and smooth- 
ing. This technique can be used with more 
limited amounts of experimental data. 

1 Hydraulic engineer, Southeast Watershed Research 
Laboratory, Agricultural Research Service, U.R. Depart 
mentof Agriculture, Athens, Ga. 30601. 


ved by extension to continuity of the 
are then transformed to smoothing 


This publication draws together and expands 
the author’s previous work on sliding polyno- 
mials. The construction of sliding polynomials is 
reviewed and analyzed. Common elements and 
concepts can be found in a hierarchy of polyno- 
mials. Tracing these common elements through 
three levels of complexity will be useful to 
researchers in many disciplines who must pro- 
cess observational or experimental data. Realiza- 
tion of the full potential of the sliding polynomials 
will aid in choosing interpolation and smoothing - 
techniques best suited to specific data processing - 
requirements. 

Understanding the hierarchial structure of the 
polynomials can begin with review of simple 
linear interpolation. Extension of this simple 
form to piecewise linear functions establishes a 
basis for understanding more complex forms. 
Addition of a requirement for mathematical 
continuity of a first derivative is a logical next 
step. Further extension follows with require- 
ment for mathematical continuity of a second 
derivative. All of these three piecewise forms are 
utilitarian. While a basic knowledge of mathe- 
matical and numerical analysis is assi 1 "’"' ) 



phasis is placed on explanation to aid the reader 
in selection of the form suited to his purposes. 

INTERPOLATION 
WITH POLYNOMIALS 

Interpolation is the estimation of values be- 
tween given discrete values of a property. Some 
function must be assumed. This function must 
give values of the property dependent upon 
values and locations of the given discrete points. 
A one-to-one correspondence exists between 
mathematical degrees of freedom in the function 
and the number of given discrete points. Data 
points through which the function is required to 
pass will be called base points. 

In the following discussion, basic polynomial 
properties and the functional forms of sliding 
polynomials are explained. Then, the use of these 
concepts in interpolation is described. 

Basic Polynomial Properties 

Polynomials of any degree could be assumed 
for interpolation. Beyond the fifth power, how- 
ever, the algebraic manipulations become tedious 
and error prone. A fifth-power polynomial, 
evaluated from six or more coordinate points, is 
sufficient to produce smoothly continuous rates 
of change through experimental data. Therefore, 
this presentation does not go beyond fifth-power 
polynomials. 

The fifth-power polynomial has the conven- 
tional form 

y~a+bx+cx: 1 +dx''>+ex*'\-fxf> ( 1 ) 

Coefficients a through / represent six mathe- 
matical degrees of freedom which can be satisfied 
by reference to six given discrete points. These 
concepts are developed under "Six-Point Inter- 
polation.” 

The first derivative of equation 1, the rate of 
change of y, is 

dy/dx=b+2cx+ , iilx i +4ex <+5/r ( 2 ) 

Higher order derivatives are given in equa- 
tions 3 through 6. 


d 2 y/dx' 2 =2c+Gdx+\2ex + 20 /x '1 ( 3 ) 

d' , y/dx 3 —6d-{-24ex-\-Q0fx‘ i (41 

rf 4 y/rfr , =24e-H20/i. ( 5 ) 

tfV<te 5 =120/. (6) 


Note the fundamental continuity of form in 
equations 1 through 5. Any one of these can be 
used as an interpolation function. Properties of 
the function for estimating rates of change are 
given by derivatives of order higher than the 
chosen functional form. Properties of the func- 
tion for integration are given by derivatives of 
order lower than the chosen functional form. For 
example, equation 5, the straight line, has a 
constant rate of change given by equation 6 and a 
parabolic rate of accumulation given by equation 
4. In the discussions of two-point, four-point, and 
six-point interpolation, equations 5, 3, and 1 are 
used to develop special polynomials. They act as 
flexible splines that pass through extended num- 
bers of coordinate points. Since symmetry is 
required, only polynomials with even degrees of 
freedom can be used in this method. 


Sliding Polynomials 

Different terminology may be used to describe 
the mathematical techniques of passing a contin- 
uously transhaping function from left to right 
through some set of discrete values, or base 
points. As the function moves to the right, base 
points must drop out on the left side of the 
function span, and new points must be picked up 
on the right side of the function span. This 
concept can be called continuous transforming, 
flexible spline, or sliding polynomial. The term 
‘‘sliding polynomial” is considered most suitable 
in the development that follows. 

First, consider the schematics in figure 1. They 
are sufficient for a simple situation where inter- 
polation is always within one interval between 
base points. 

Note the set of coordinate points FBO through 
FJ36 in figure la. These coordinate points, spaced 
at uniform intervals, will be used as base points 
through which interpolation polynomials are re- 
quired to pass. The designation YBO, for exam- 
ple, means a base point in y at a value of x equal 
to zero. Uniform spacing of the base points is not 
absolutely necessary, but it does greatly simplify 
the algebraic forms. 

Now observe in detail the interval between 
YB2 and YB3. Figure lb is the schematic 
representation for passing a straight line be- 
tween the two points. Interpolation can be made 
on this straight line. Figure lc shows a cubic 
equation passing through FBI, FB2, FB3, and 
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a DEFINITION OF BASE POINTS 


VD2 YB3 

O O 

y- a + bx 

b. TWO- POINT INTERPOLATION 

Y02 YB3 

Y0I „ - - ** ^ ^ ^ 

^ v xYB4 

y= a + bx + cx 2 + dx^* 

c. FOUR“ POINT INTERPOLATION 



y = a + bx + cx 2 + dx 3 4* ex'* + fx 5 

d. SIX-POINT INTERPOLATION 

FIGURE 1 —Definition of interpolation systems. 

754. Interpolation is to be limited to the interval 
between 752, and 753, the same interval cov- 
ered by the straight-line example. Figure Id 
shows a fifth-power equation passing through six 
points. Interpolation again is to be limited to the 
interval between YB2 and 753. 

A more realistic, and more complicated, prob- 
lem is interpolating in adjacent intervals (fig. 2). 
Figure 2a shows interpolation on two adjacent 
linear segments. Visualize a point of interpola- 
tion on the x-scale moving from YB2 to 753. 
Values of y can be read from the straight line 
connecting YB2 and 753. When the point passes 
through 753, a transformation is necessary. 
Point YB2 is dropped on the left, and point YBA is 
added on the right. As the x-poinl moves through 
x=3, interpolation changes abruptly from the 
line through 752 and 753 to the line through 
753 and YBA. 

A smoother transition through 753, utilizing 
an arc through 752 and 753 and one through 
753 and 754, is shown in figure 2b. The 
technique to suppress the angle between linear 
segments at 753 is developed under “Four-Point 
Interpolation.” It is based on dropping out 751 


on the left and picking up 755 on the right. The 
phasing out and phasing in of the base points is 
gradual, instead of abrupt, as in the linear 
segments. 

Figure 2c is an extension of concepts in figure 
2b. The higher order polynomial phases out 750 
and phases in 756 as the x-point leaves the 

752- 753 arc and passes through x=3 into the 

753- 754 arc. Special forms for this method are 
shown under "Six-Point Interpolation." 

Two-Point Interpolation 

Two-point interpolation is the familiar and 
conventional linear interpolation. The method is 
intuitively understood and used by many per- 
sons. This intuitive understanding is used here as 
the basis for presentation of more complex 
polynominal methods. To prepare for this pres- 
entation, the familiar linear interpolation is cast 
into algebraic forms which will be needed later. 
Only two base points at a time are used in linear 
interpolation. In figure la, the two points are 
752 and 753. The two parameters in the 
associated equation can be evaluated in terms of 
the base points. The matching of number of 
parameters and number of base points should be 
noted. 



a. CONNECTED LINEAR SEGMENTS 



b. SLIDING PARABOLA 



c, SLIDING FOURTH POWER POLYNC 

FIGURE 2. — Definitions for connected arcs. 
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Define a new scale variable z, with z=0 at 
x—2 t and z=l at x=3, An equation to define the 
assumed linear interpolation function between 
752 and 753 can now be written: 

y — a+bz, (7) 

At z=0, y-YB2, and z-l, */=753. Substi- 
tuting these coordinate values of the two base 
points in equation 7 yields the simultaneous 
equations 

YB2 =a (8) 

YB3~a+b (9) 

From equations 8 and 9, a=752, and 6 = 753 
-752. Equation 7 may now be written in terms of 
the base points as 

y=YB2.MYBb-Y&)z=(\-z)Ym+*YBZ. ( 10 ) 

The solution for equations 8 and 9 may also be 
written in the form of the matrix-vector equation 


1 0 


YB2 


a 

-1 i 


YB2> 


b 


So-called interpolation coefficients can be de- 
veloped from special values of 752 and 753. 
First, set 752 equal to 1 and 753 equal to zero. 
Expansion of equation 11 gives a = l and 6=-l. 
The special form of equation 10 results in 

1“*. U2) 

Since 752 was set to unity, y a can be called the 
unit contribution of a base point at z=0, the left 
one of two base points. Next, set 752 equal to 
zero and 753 equal to 1. Expansion of equation 
11 gives a=0 and 6 = 1. The special form of 
equation 10 now is 

(13) 

Now, y b can be called the unit contribution of a 
base point at z= 1, the right one of the two base 
points. 

Equations 12 and 13 can be evaluated for 
values of z as in table 1. The values of y a and y b 
are normally designated as CO and Cl, meaning 
the interpolating coefficient to be multiplied by 
the base value at z^O and by the base value at 

— . i rp l „ _ , , 

meaning is expressed on the 
n F or example, two-point 
-fourth of the i-dis- 
ritten directly, with 

- - — x, ao 

2f=l0.75)m-H0.25)ra3. (H) 

Equation 14, of course, is the familiar and 


conventional form for linear interpolation. 

If the two adjacent linear segments in figure 2a 
are considered, equations can be written directly 
for interpolating y , on the left segment, and y , on 
the right segment: 

)/ =YB2+I YB3- YI)2)z, ( I S) 

and >/ =YBll+[YM-Ym)z. (Hi) 

With z— 1, equation 15 gives ?/, = 753. Also, 
with 2=0, equation 16 gives y ,— YB'S. This 
simply confirms that both segments pass through 
753, and the interpolation is continuous in the 
sense that there is no gap at the base point 753. 
However, dy , ldz-YB3-YB3 and <iy>/dz— YB4 
-753 show that the slopes of the segments are 
different, constant, and not continuous. 


Table 

1 — Two -point interpolation 

coefficients 

Z 

Vn 

V/j 

0 10 

0 90 

0.10 

.25 

.75 

.25 

.50 

.50 

.50 

75 

.25 

.75 

90 

.10 

90 


Four-Point Interpolation 

In geometrical concept, two-point interpola- 
tion can be visualized as a point being linearly 
displaced in (/-scale as it moves along the 2 -scale 
from one base point to the next. This two- 
directional displacement generates a locus of 
interpolated values which lie along a straight 
line. Note particularly the concept of a point in 
continuous displacement generating a straight- 
line interpolation locus, and note that the line is 
one mathematical order up from the point. 

In tour-point interpolation, the geometrical 
concept is a function that will be linearly trail- 
shaped in y-scale as a point of interpolation 
moves along the 2 -scale. This function is up two 
orders from the point, and is, therefore, a 
parabola. The locus of interpolated values is one 
Older higher than the function being transhaped 
and is a cubic polynomial. 

1 he geometrical concept on which four-point 
interpolation is based is further developed as 
follows: Visualize a parabola that passes through 
base points 751, 752, and 753 in figure 1. Next, 
visualize the parabola through 752, 753, and 
754. As the point of interpolation moves along 
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the x-scale from YB2 to YB3, the interpolated 
value is to be read from a flexible parabola that 
transhapes linearly from the left parabola to the 
right parabola. 

The analytical concept of four-point interpola- 
tion requires specification of four mathematical 
constraints on which solution for the four coeffi- 
cients of the cubic polynomial, equation 3, are to 
be based. The mathematical specifications are 
based on commonality. The cubic equation is 
required to be common to the parabolas at the 
ends of the interpolation arc. This means simply 
that both the parabolas and the cubic pass 
through YB2 and YB3 (fig. lc). In addition, the 
parabolas and the cubic are required to have the 
same slope, or the same first derivative, at YB2 
and YB3. 

Next, specify a scale variate z, which takes 
values 0, 1, 2, and 3 from left to right in any 
four-point set of base points. Define the left 
parabola by 

, y—u+vz+wz- (17) 

and d 1 i[/dz=v+2wz. (18) 

Solve for coefficients u, v, and w by passing 
through the three points 2=0, l y=YB\; 2=1, 
{ y—YB2; and z— 2, t y~YB3, giving the values 



m =ra 1, 

(19) 


v =-<3/2)17*1+ 2 YB2-tt/2)YB\l, 

(20) 

and 

w = ( 1 n )YBU YB2 +(1/2) YM 

(21) 

Next define the right parabola by 



1 >l/~u+vz+wz 

(22) 

and 

d iy/dz—v+2wz 

(23) 


Solve for new coefficients u, v, and w by 
passing the parabola through points 2 = 1, 
■>y-YB2; 2=2, ,y = YB3; and 2=3, ,y~YB4. 

The solution is 


u— 3' KB2-3- FJ33+ KB4, (24) 

v=-(5/2)yU2+4KB3-(3/2)yi?4, (25) 

and ?c=(l/2)yB2-yfl3+(l/2)yS4. (26) 

Define the cubic equation by 

y=a+bz+cz -bdr 1 (27) 

and dy/dz=zb+2cz-\-3dz 2 (28) 


The four mathematical constraints expressing 
commonality between parabolas and cubic can 
now be stated as four simultaneous equations 
(29-32). Equations 17 and 18, with coefficients in 
equations 19, 20, and 21, are set equal to 


equations 27 and 28 for 2=1. Equations 22 and 
23, with coefficients in equations 24, 25, and 26, 
are set equal to equations 27 and 28 for 2=2. The 
resulting four simultaneous equations are 


YB2=a+b+c+d, (29) 

[\/2)YB3-{\/2)YB\=b+2c+3d, (30) 

YD3=(i+2b+4c+-B<l, (31) 

ami (l/2)yit4-(l/2>y/12=6+4c+12</. (32) 


The solution to these equations is 

a =2 • y/n-3- KZ12 + 3- YB3- YM . (33) 

(,=-4.y2n-H19/2)KQ2-.8-yB3H-(5/2)yS4, (34) 

c— (5/2) YB l-l ■ yfl2-|-(13/2) YB3-2YM, (35) 


.mil d=-(l/2)yfll-f(3/2>yB2-(3/2)yS3-f(l/2)yfl4. (36) 

This solution can be expressed concisely in the 
matrix-vector equation 


4 

- 6 

6 

- 2 


YB1 


« 

- 8 

19 

-16 

5 


YB2 


b 

5 

-14 

13 

- 4, 


YBS , 

— 

c 

- 1 

3 

- 3 

1 


YM 


d 


Next, consider the property of the locus of 
interpolated values at the common point of two 
arcs. Specifically, in figure 2b, arcs YB2-YB3 and 
YBH-YB4 obviously join at YB3, Now at YB3, 
interpolation arc YB2-YB3 ends on the parabola 
through YB2, YB3, and YB\. However, interpo- 
lation arc YB3-YBA takes off from this same 
parabola. The four mathematical specifications 
require the interpolation cubics to be common 
and tangent to the parabolas at the ends of the 
arcs. The cubics of the YB1-YB4 set and the 
YB2-YB5 set are, therefore, both common and 
tangent to the same parabola at YB3. The cubics 
of adjacent interpolation arcs are thus common 
and tangent to each other. 

Thus, it has been demonstrated that four-point 
interpolation based on a sliding, transhaping 
parabola generates values on a cubic polynomial. 
Adjacent cubic arcs are continuous and have 
continuous first derivatives as the interpolation 
locus passes through the base points. 

Figure 3 is an example of four-point interpola- 
tion. (Computer program A-l, appendix A, gen- 
erated the interpolated values.) Remember, 
when looking at figure 3, that each interpolation 
arc is a different cubic polynomial, with the 
ensemble of cubics having the unique property 
that they are common and tangent at the base 
points. 
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FIGURE 3 — Four-pomi interpolation and integration. 

In summary, four-point interpolation estab- 
lishes a form-free ensemble with two orders of 
continuity. Generally, the interpolation locus will 
pass through the base points without undesirable 
oscillation. In some situations, with an abrupt 
change in slope of the base point data, a tendency 
to oscillate may develop. 

Four-point interpolation can be reduced to 
coefficients that operate directly on the base 
points, just as the two-point system was re- 
duced. Consider the special case of 751=1 and 
752 = 753 = 754 =0. Then, from equations 33-36 
or equation 37, a= 2, f>=-4, c= 5/2, and d= 
"1/2, and the interpolation cubic is 

&„=2-4^-H5/2)i jl -(l/2)2>. (38) 

Equation 38 can be evaluated for selected 
values of z, producing a coefficient that gives the 
unit contribution of the base point YB1 to the 
interpolated value. Since YB1 is at a z-value of 
zero, call this coefficient CO. 

^Now use the special case 752=1 and YB1 
— 753=754=0. The interpolation cubic from 
equations 33-36 is 

3/ 1 = 3+ ( 19/2)2-72- + 1 3/2)z' . (39) 

Evaluate this for the same selected z-values 
used with equation 38. Call this set of coefficients 


Cl, since 752 is at a z-value of 1. 

Similarly, evaluate C2 and C3. An interpolated 
value is the summation of the products of Lhe unit 
operators with any real sot of 751 through 754 . 
This product summation is 

y=y {) +y i +y„+y i-CO-YBl+Cl 752 + 62 753 

+53-754 (40) 

Table A-l (appendix A) is a listing of the four 
interpolating coefficients for each one-hundredth 
of the interpolating arc from z=l to z~ 2. These 
coefficients were generated by program A-2. 

The area under any function is given by the 
integral of the function across a range of x (or z). 
The integral of the interpolation cubic of the 
four-point method is given in equation 41. Inte- 
gration, just as interpolation, must be limited to 
the central interval of the four-point set: 

^ ydz—az+(b/2)z'- ■\-{c/ , i)z '+(d/A)z'-\-C=A (41) 

The constant of integration, C, can be evalu- 
ated by noting that the area under the curve, A, 
is zero when z=l; that is, at the left end of the 
interpolation arc. For this value of z, and with 
equations 33-36, C is given by 

C- ( l/24)( - 17- 751 + 5- 7712- 1 9 • YDS + 7 • 754 ) - (42) 

The area under any segment of the interpola- 
tion arc can now be calculated with equations 
33-36, 41, and 42. It is also possible to reduce 
four-point integration to a coefficient-multiplier 
scheme analagous to equation 40 for interpola- 
tions. lable A-2 is a listing of such coefficients, 
also generated by program A-2 (appendix A). 
Now, note that four-point integration is also 
included in program A-l, The area under the arcs 
of the interpolation cubic are plotted in figure 3. 
Remember that each arc starts with zero area, 
lhe total area to the left of any point is thus the 
sum of the partial area in the arc plus all full arcs 
to the left. 

The coefficients in table A-2 are based on a 
scale variate z with the base points at unit 
intervals of z. If the base points are at a scale s in 
the original data, given by . then 

ds—kdz . An integral ydz can be changed to an 
integral yds—kydz by simple change of variable, 
This change of variable is included in program 
A-l, but must be performed as a separate 
operation following use of table A-2. 
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Six-Point Interpolation 


The development of the method of six-point 
interpolation has been given in detail previously 
(Snyder 1967a). The geometrical and analytical 
concepts are analagous to those used in develop- 
ment of the four-point method in the previous 
section. 

Note briefly that six-point interpolation re- 
quires functions two orders higher than four- 
point interpolation, which in turn was two orders 
higher than the functions of two-point interpola- 
tion. Specifically, six-point interpolation is based 
on a sliding, transhaping, fourth-power polyno- 
mial, two orders up from the sliding parabola. 
Each arc of the interpolation locus is given by a 
fifth-power polynomial. 

Six mathematical restraints are required. In 
addition to the four used for four-point interpola- 
tion, commonality of function and first deriva- 
tive, commonality of the second derivative with 
' the second derivative of the linearly transhaping, 
fourth-power polynomial is required at the ends 
of the interpolation arcs. The junctures of the 
arcs at the base points are mathematically 
continuous and have continuous first and second 
derivatives. 

Details of the solution of the coefficients of the 
six-point method are not repeated here. The 
solution for the coefficients in equation 1 is 
conveniently expressed by the matrix-vector 
equation 


(1/24) 


432 -2040 4080 -4080 2040 -408 


YB 0 


a 

-918 4436 -8752 8712 -4346 868 


YB\ 


b 

765 -3754 7414 -7366 3661 -730 


YB2 


c 

-313 1551 -3078 3058 -1521 303 


YBS 


d 

63 -314 626 -624 311 -62 

i 

YB\ 


e 

-5 25 -50 50 -25 5 


YBS 


/ 


H3) 


In the previous section on four-point interpola- 
tion, it was pointed out that the method could be 
extended to integration of experimental data. An 
advantage of the six-point system is that it can be 
converted to differentiation, leading to deter- 
mination of smooth rates of change of experi- 
mental data. The rate of change of values in 
equation 1 is given directly by equation 2. With 
coefficients b through / given by equation 43, 
rates of change can be computed by direct 
substitution in equation 2. An advantage of this 
method of interpolative differentiation is that the 
first derivatives are smoothly continuous be- 
cause the second derivatives are continuous. 


F'igure 4 shows six-point interpolation through 
and the calculated first derivative of the same 
data as in figure 3. The interpolative values and 
derivatives were computed with program A-3 
(appendix A). 

Six-point interpolation and differentiation can 
be reduced Lo multiplying coefficients in a man- 
ner similar to interpolation and integration with 
coefficients as presented in the four-point 
method. Interpolation coefficients for each one- 
hundredth of the interval of the interpolation arc 
from z = 2 to z=3 are given in table A-3 (appen- 
dix A). Those coefficients were generated by 
program A-4. Table A-4 gives coefficients for use 
in differentiation. Calculation of these coeffi- 
cients is included in program A-4. 

It should he noted that slopes given by the 
differentiation coefficients of table A-4 are in 
terms of the scale variate z. Given that the scale 
of original data, s, is related to z by s=kz+z„, 
then ds—kdz. A slope dy/dz can be changed to a 
slope of original data by simple change of 
variable, yielding dy/ds~{\/k) dy/dz. This scale 
conversion is included in program A-3, but must 



FIGURE 4. — Six-point in lorpo lotion and differ 




hi 1 performed as a separate operation following 
use of Lable A-4. 

SMOOTHING OF 
EXPERIMENTAL DATA 

Smoothing of experimental data is essentially 
the application of the method of least squares to 
effect a transformation from recorded data to a 
set of base points for interpolation. Data may he 
recorded at nonuniform intervals. Transforma- 
tion produces base points at uniform intervals. 
Data may include mensuration errors. These 
errors may be smoothed by averaging during the 
transformation. Successful averaging requires 
that base points be fewer in number than data 
points and also requires a configuration of data 
points conducive to production of good averages. 

The location and spacing of the base points is a 
problem in design. No hard rules can be stated. 
The heart of the problem is the amount of 
averaging that is desired. With a goodly amount 
of data in a compact mass, a high degree of 
averaging can be used. The base points would 
likely be one-fifth to one-tenth of the number of 
data points. With limited data, strung out in 
range, the base points would likely be about 
one-half the number of data points, resulting in a 
limited degree of smoothing. 

In the following examples, a high degree of 
smoothing will be evident. The placement of 
points is usually dictated by the particular 
problem. Abase point would usually be placed at 
an interface or boundary, such as the surface of 
ground or water. In time series one point would 
usually be placed at time zero. Spacing of the 
remaining base points through the data range is 
obviously dictated by the degree of averaging 
and hence the number of points. 

Two-Point Smoothing 

Two-pomt smoothing is identical to evaluation 
of connected linear segments (Snyder 1967b), 
Equation 14 can be generalized as 

y=CQ-YBQ+CbYB\. (44) 

In equation 44, the meaning is limited to one 
interpolation interval. If several adjacent inter- 
vals are considered, equations 45-47 result: 


y, -CO-FBO+CM’Bl 

(45) 

V 2~ 

C0-YB1 + CI-YB2 

(46) 

y^~ 

C0-YB2+C1YB3 

(47) 


In equal ion 45, y is any value lying on the line 
from YB 0 to Yll 1, In equation 46, //, is any value 
lying on the line Irom YD l to YB2, and so on. The 
interpolation coefficients CO and (A are known 
functions of scale variate z, as in table 1. 

Assume that experimental data consist of 
simultaneous coordinate values ol .r and y lor a 
series of observations. Thou, consider a series of 
base points to be substituted for the data points. 
It is necessary to specify the location of the base 
points within the .v-seale of the data. Normally, 
one has some idea of the number of base points 
desired to span I he data in the .r-scalo. Setting 
the location of tin* base points in .e-scale automat- 
ically sets a ./-scale value for each data point since 
the base points are always spaced at unit interval 
in z. 

Referring now to equations 45-47, note that 
they are linear in form. Note also thal given 
experimental values of y and given a computed 
z-variate for each observation, there is an equiv- 
alent of a series of observational equations as 
used in multiple regression. CO and Cl take the 
position of x - ,, and a: , in regression, and YU 0 and 
YB 1 lake the posit ion of regression coefficients, 
conventionally designed and b |. There is no 
regression intercept normally designated a. It is 
a simple matter, with conventional least -squares 
methods or programs, to compute YBO and YB 1 
as regression coefficients. A transformation from 
experimental data in x-y coordinates to base 
point data, in z-YB\ coordinates, is thus accom- 
plished. Note especially that in the configuration 
of equations 45-47, VB1 and YB2 must take 
values on two adjacent intervals. The continuity 
of Lhe smoothing operation is thus satisfied with 
no gaps at the base points. 

Table U-l (appendix B) shows the details for 
organization of 30 data points into format for 
regression transformation. The data poinl.s and 
the result of least-squares transformation in 
two-point smoothing are plotted in figure 5. 
While the location of base points at .c=0.4, 0.6, 
0.8, and 1.0 is arbitrary, the values of these base 
points are wholly objective and the result of 
least-squares fitting. Note specifically in figure 5 
the overall fit of the linear segments to the rough 
crescent shape of the experimental data. Since 
the two internal base points are each required to 
be common to two segments, the fit of a single 
segment to the data points in that segment is 
secondary. 
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FIGURE 5 - — Smoothing by two-point method. 

Four-Point Smoothing 

Four-point smoothing is a simple extension of 
two-point smoothing. Equation 40, modified for 
adjoining arcs, is extended in equations 48-151, 
analagous to the two-point equations, 45-47. 

H ] ~C0-YB0+Cl-Ym + C2-YB2+C3 YB3 [-181 

C0>YB2+CbYBZ+C2-YB3+Ca-YlM (49) 

U.= CQ-YB2-\-C\-YB3 

+C2-YB4+C3-YB5 (50) 

//.,= CQ-YIM+CX-YM 

-i-C2-YB5+Gi-YB(t (51) 


Values for the coefficients are given in table 
A-l. Remember, in equations 45-47 each base 
point is linked to two adjacent intervals, YR1 and 
YB2 in those equations are each linked by CO and 
Cl. Note in equations 48-51 that an internal base 
point such as YB3 is linked to four interpolation 
arcs. The least-squares value of YB3 is depend- 
ent on experimental data in four intervals. 
Exterior base points are dependent on 1, 2, or 3 
intervals. 

Table B-2 (appendix B) shows organization of 
experimental data in format for four-point 
smoothing. The same data are used as for the 
two-point smoothing example, except the a-scale 
is shifted to the left to put the first base point at 
X'=0.4. The z-scale is not affected. Program B-l 
(appendix B) is used to perform the least-squares 
transformation. This program also generates the 
locus of interpolated values which runs through 
the derived base points and generates the inter- 
grals of the interpolation arcs. 

Figure 6 shows the results of the four-point 
smoothing. An interpolation locus of smoothly 
joined arcs is generated that runs continuously 
through the data. Comparison of figures 5 and 6 



FIGURE 6 — Smoothing by four-point rnefhod. 


shows that the two-point and four-point methods 
produce base points not greatly different. The 
main difference is the creation of arcs by the 
four-point method, smoothly joined with angles 
eliminated. Values of the derived base points can 
be compared in table 2. 

The integral of the interpolation arcs is also 
plotted in figure 6. This integral can be thought 
of as the smoothed summation of the experi- 
mental data. 

Six-Point Smoothing 

Six-point smoothing is a simple extension of 
methodology already presented. Equations 52-55 
are analogous to equations 45-47 and 48-51. 

1/ , = CO • YBO +CI-YIH+C2' YB2 

+ C3 • rail + C4 ■ YB4 + C5 • YB5 (52) 

)/. =C0- Ym+Cl- YB2+C2- YB'A 

+C3-YB4+Ci-YBr)+C5-YB6 (63) 

H =CQ-YB2+Cl-YB3+C2-YB4 

+C3'KB5-)-C4 YBGYCS'YBl (54) 

n=C0-YB\\+C\-YB4 + n YB§ 

+ C3 • YBti -r- CA ■ YB7 + C5- YD8 (56) 

In equations 52-55, note that an interior base 
point will depend upon values of experimental 


TABLE 2. — Comparison of derived base points 


Two Point 

Foiir-Pomt 


.Six-Point 

X 

Ylh 

X 

rat 

X 

YBi 





0.4 

-4.14772 



0.4 

-1,72742 

,6 

-1.55935 

0 4 

1.52418 

.6 

1,42889 

.8 

1.417502 

.6 

8.75141 

.8 

3.59534 

1.0 

3.69252 

.8 

8.81551 

1.0 

3.67674 

1.2 

3,68015 

1.2 

2.57954 

1.2 

2.44932 

1.4 

2.45502 



1.4 

- .42235 

1.6 

-2.29737 





1.8 

-20.17608 
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data in six adjacent intervals by linkage through 
the coefficients. Exterior base points may de- 
pend on five or fewer intervals of data, 

Table B~3 (appendix B) shows organization of 
experimental data in format for six-point smooth- 
ing. Program B-2 is used to perform the least- 
squares transformation. Following transforma- 
tion of the experimental data to a smoothed set of 
base points, program B-2 computes a locus of 
interpolated values and a locus of differentials of 
the interpolation ares. Figure 7 is a plot of the 
experimental data and derived interpolation and 
differentiation loci. 



Figure 7 contains the same experimental data 
as figures 5 and 6. It can be seen that the 
interpolation arcs of the four-point and six-point 
methods are almost identical. Numerical com- 
parison of the derived base values is shown in 
table 2. 

Figure 7 also shows the differentiation of the 
interpolation arcs. This locus of values can be 
regarded as the smoothed rate of change of the 
experimental data. Note the smoothness of this 
locus through the base points and the near 
constant value of slope from about £=1.14 to 
£= 1 . 2 . 

ATTRIBUTES OF 
SLIDING POLYNOMIALS 

Figures 5, 6, and 7 demonstrate that various 
interpolation methods can be used to smooth 
experimental data. Smoothed values, whether 
read from two-point, four-point, or six-point 


interpolation loci, are not extremely different. 
Other methods would also produce acceptable 
interpolation loci. Kimball's (1974) method would 
require a means to separate the spectral portions 
of noise and information in the data. 

The identifiable advantages of sliding polyno- 
mials are in the secondary characteristics of their 
loci. The four-point and six-point methods are 
smooth because their derivatives are continuous 
at the juncture points of the interpolation arcs. 
The four-point method, with a continuous first 
derivative, produces an extremely smooth locus 
of integral values, since integration itself is a 
smoothing process. The six-point method, with 
continuous first and second derivatives, can gen- 
erate a smooth rate of change in experimental 
data. DuGhateau et al. (1972) also used six points 
in evaluation. In their method, however, para- 
bolic splines are required to be continuous only in 
the first derivative. Therefore, these splines 
cannot be expected to develop the smoothly 
continuous rates of change of experimental data. 
Such continuous smoothness is dependent upon a 
continuous second derivative. 

Sliding polynomials and their integrals and 
derivatives can be derived from experimental 
data by conventional least-squares techniques. 
The interpolation coefficients provide exact error 
weights for each experimental data point, re- 
gardless of the interpolation interval in which the 
data point is located. The structure of equations 
52-55, for example, shows that the evaluation 
procedure is similar to the use of smoothing 
convolutes as in Savitzky and Golay (1964). 
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appendix a.— coefficients and computer programs for 

INTERPOLATION 

TABLE A-l .— Coefficients for four-point interpolation 


1.00 

0.00000000 

1.00000000 

.uooooooo 

0 . 00 000000 

1 . 01 
1.02 

-.00490050 

-.00960400 

.99975450 

.99901200 

. 00519850 
.0 1078800 

00004950 

-.00019600 

1.03 

1 . 01 * 

-.01411350 

-.01843200 

.99779050 
.9960960 0 

.01675950 

.02310400 

-.00043650 

-.00076800 


1 . 05 
1.06 

-.02256250 

-.02650800 

.99393750 

.99132400 

,02981250 
, 03687600 

00 118750 
-. 00 169200 

1.07 

1.06 

-.03027150 

-.03385600 

. 98826450 
. 98476 BU 0 

.04428550 
.0520 3200 

-. 00 227850 
00 294400 

1 . 09 

-.03726450 

.98084350 

.06010650 

-. 00 368550 

1.10 

- .04050000 

.97650000 

, 0 685 u 00 0 

-. 00 450000 

1.11 

1.12 

-.04356550 

-.04646400 

.97174650 

.96659200 

. 07720350 
.0862080 0 

-. 0 0 538450 
-, 00633600 

1.13 

1.14 

-.04919850 

-.05177200 

.96104550 

.95511600 

.09550450 
.10 50 840 0 

-. 00 735150 
- . 00 842800 


1 . 15 
1. 16 

-.05418750 

-.05644800 

. 94881250 
.94214400 

.11493750 

.12505600 

-. 00956250 
-. 01075200 

1.17 

1.18 

-.05855650 

-.06051600 

.93511950 

.92774800 

,13543050 
.14605 200 

-.01199350 
-. 01328400 

1. 19 

-.06232950 

, 92003850 

,15691150 

-. 01462050 

1.20 

-.06400000 

.91200000 

.16800000 

-.01600000 

1 . 21 
1.22 

-.06553050 

-.06692400 

. 90364150 
. 89497200 

.17930850 

,19002800 

-. 01741950 
-. 01887600 

1 . 23 
1.24 

-.06818350 

-.06931200 

. 88600050 
. 87673600 

.20254950 

.21446400 

-. 02036650 
-. 02188800 


1.25 

1.26 

-.07031250 

-.07118800 

. 86718750 
.85736400 

.22656250 
. 23883600 

02343750 
-. 02501200 

1.27 

1.28 

-.07194150 

-.07257600 

. 84727450 
. 83692800 

. 25127550 
,26387200 

- . 02660850 
- . 02822400 

1.29 

-.07309450 

.82633350 

.27661650 

-. 02985550 

1.30 

-.07350000 

. 81550000 

.2895000 0 

- ,03150000 

1.31 

1.32 

-.07379550 

-.07398400 

• 804 4 3650 
. 793152 G 0 

. 3 0251350 
. 31564800 

-.03315450 
-. 03481600 

1.33 

1.34 

-.07406850 

-.07405200 

.78165550 

.76995600 

. 32889450 
,34224400 

03648150 
-. 03814800 



TABLE A-l .—Coefficients for four-point interpolation —Continued 


i , 35 
1 . 36 

- . 07393750 
07372800 

.75806250 

.74598400 

. 35568750 
. 3692160 0 

03981250 

-.04147200 

1 , 37 
1.36 

“.07342650 

“.07303600 

. 73372950 " 

. 7213080 0 

.38282050 
. 3964920 0 

-.04312350 

-.04476400 

1.39 

“.07255950 

.70872850 

.41022150 

04639050 

1.40 

“.07200000 

.69600000 

.42400000 

-.04800000 

1.41 
1 . 42 

- , 07136050 
“.07064400 

.68313150 

►67013200 

.43781850 
.4516680 0 

-.04958950 

05115600 

1.43 

1.44 

“ .06985350 
“.06899200 

.65701050 

.64377600 

.46553950 
. 479^240 0 

-. 05269650 
-.05420800 

1.45 

1.46 

1.47 
1 . 48 

- .06806250 
“ .06706800 
- . 06601150 
“ . 06489600 

.63043750 

.61700400 

.60348450 

.56988800 

.49331250 
. 50 71960 0 
.52106550 
. 5349120 0 

- . 05568750 
-.05713200 

- . 05 853850 
-. 05990400 

1.49 

- . 06372450 

.57622350 

•54872650 

-.06122550 

1.50 

0 6250 300 

.56250000 

.56250000 

-.06250000 

1.51 

1.52 

“.06122550 

-.05990400 

.54872650 
. 53491200 

.57622350 

.58988800 

-.06372450 

-.06489600 

1.53 

1.54 

-.05853850 

-.05713200 

.52106550 

.50719600 

.60348450 

.61700400 

06601150 

-.06706800 


1.55 

1. 56 

- . 0 5568750 
-.05420800 

.49331250 
.4794240 0 

.63043750 
.6437760 Q 

06806250 

-.06899200 

1.57 

1.56 

-.05269650 
- .05115600 

.4655395 0 
.4516680 0 

. 65 70105 0 
.67013200 

-.06985350 

-.07064400 

1.59 

- .04958950 

,43781850 

. 6831315 0 

-.07136050 

1.60 

- .0 48 00000 

*42400000 

. 6960 u 00 0 

-.07200000 

1.61 

1.62 

-.04639050 

-.04476400 

.41022150 
. 39649200 

.70872850 

.72130800 

-.07255950 

-.07303600 

1.63 

1.64 

-.04312150 

-.04147200 

. 38282050 
.36921600 

.73372950 

.74598400 

07 342650 
-.07372800 


1.65 

1.66 

-.03981250 
-. 0 3814800 

.35568750 
. 34224400 

. 75806250 
.76995600 

07393750 

-.07405200 

1.67 

1.68 

- . 0 3648150 " 
-.03481600 

.32889450 
. 31564800 

.78165550 
.7931520 0 

-.07406850 

07398400 

1.69 

-.03315450 

.30251350 

.80443650 

-.07379550 



TABLE A-l. — Coefficients for four-point interpolation— Continued 


1.70 

- . 03150000 

. 28950000 

. 8155 0000 

07350000 

1.71 

1.72 

-.02905550 

-.02822400 

. 2766165 0 
. 26387200 

. 82633350 
. 83692800 

-.07309450 

-.07257600 

1.73 
1. 7 4 

- .02660850 
- . 0 2501200 

.25127550 
. 2388360 0 

. 84727450 
. 85736400 

07194150 

-.07118800 


1. 75 
1. 76 

- .02343750 
- .02188900 

.22656250 

.21446400 

, 86718750 
. 87673600 

-. 070 31250 
06931200 

1.77 

1.78 

- .02036650 
- .01887600 

. 20254950 
. 19082800 

,88600050 
. 89497200 

-.06818350 

-.06692400 

1.79 

- .01741950 

. 17930 850 

.90 36415 0 

-.06553050 

1.80 

- .01600000 

. 16800000 

. 912G 0 00 0 

-.06400000 

1.81 

1.82 

-.01462050 

-.01328400 

.15691150 
. 14605200 

. 92 00 3850 
. 92774800 

-.06232950 
-. 06051600 

1. 83 
1.84 

- .01199350 
- . 0 107520 0 

.13543050 

.12505600 

. 93511950 
. 94 214400 

-.05855650 

-.05644800 


1.85 

1.86 

- .00956250 
- .00842800 

.11493750 

.10508400 

. 94881250 
. 9551160 0 

-.05418750 

-.05177200 

1.87 

1.38 

- .00735150 
- . 00633600 

.09550450 

.08620800 

.96104550 
. 96659200 

-. 04 919850 
-.04646400 

1.89 

- .00538450 

.07720350 

.97174650 

-.04356550 

1.90 

- .0045000 0 

.06850000 

.97 65 0 00 0 

-.04050000 

1.91 

1.92 

-.00368550 
- .00294400 

. 06010650 
.05203200 

.98084350 
, 98476800 

-.03726450 

-.03385600 

1.93 

1.94' 

-.00227850 
- , 0 0 1 69 2 0 0 

,04428550 
.0 3687600 

. 96826450 
.99132400 

-.03027150 

-.02650800 

1.95 

1.96 

1.97 

1.98 

-.00118750 
- .00076800 
-. 00043650 
-. 000 19b u 0 

. 02981250 
.02310400 
. 01675950 
. 01078800 

.99393750 
. 99609600 
. 99779050 
.99901200 

-.02256250 

-.01843200 

-.01411350 

-.00960400 

1.99 

- .000 04950 

.00519850 

. 99975 150 

-.00490050 

2.00 

O.OuOOOOOO 

-.00000000 

1. 00 000 00 0 

0 . OOOOOOOO 
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TAULE A-2.— Coefficients for four-point integration 


JL. 

vXQ 

Cl 

C2 

C3 



* ) • 

1.00 

-.00000000 

OODOUQuO 

0. Ou OQQGOO 

. Gu 000000 

1.01 

1.02 

“.00002467 

-.00009735 

.00999917 

.01999339 

.00002566 

.00010527 

00000017 
00 000131 

1.03 

1.04 

Q0U2161O 

-.00037899 

.02997780 

.03994763 

. 00024270 
. OC U44171 

00000440 
00 00 10 35 


1.05 

-.00058411 

.04969818 

. 00070599 

00 00 2005 

1. 06 

-.00082962 

, 05982486 

. Ob 10 3914 

OU 003438 

1.07 

-.00111 367 

.06972317 

.00144466 

00005417 

1.08 

-.00143445 

. 07958869 

.00192597 

00008021 

1.09 

-.00179020 

.08941710 

.00 248640" 

00011330 


1.10 

-.00217917 

♦ 09920417 

.00312917 

Q0015417 

1.11 

-.00259963 

, 10894574 

.00385743 

OU 020 353 

1.12 

-.00304992 

. 11863776 

.00467424 

G 0 U 2620 8 

1.13 

-.00352837 

.12827627 

.00558256 

-. 00033047 

1.14 

-.00403335 

.13785739 

.00658527 

-. 00 040931 


1.15 

1*16 

-.00456328 

-.00511659 

.14737734 
. 15683243 

.00768516 
. 00 888491 

00049922 

-.00060075 

1. 17 
1.18 

-.00569173 

-.00628722 

. 1662190*+ 
.17553366 

.01018713 

.01159434 

00 071443 
00084070 

1.19 

-.00690157 

.1847728 7 

.01310696 

-. 00 098027 

1.20 

-.00753333 

.19393333 

.01473333 

-. 00113333 

1.21 

1.22 

-.00818110 

-.00884349 

. 20301180 
.2120051 3 

.01646970 

.01832021 

QO 130040 
00 148185 

1.23 

1.24 

-.00951913 

-.01020672 

. 22G91024 
.22972416 

. 02028693 
.02237184 

00 167803 
00188928 


1.25 

1.26 

-.01090495 

-.01161255 

. 23844401 
. 24706699 

.02457682 

.02690367 

-. 00 211589 
00235811 

1. 27 
1.28 

-.01232830 
- . O'l 3050 99 

. 25559040 
« 26*+01163 

.02935410 

.03192971 

00 261620 
-. 00289035 

1.29 

-.01377943 

. 27232814 

.03463203 

00318073 

1.30 

-.01451250 

. 28053750 

.03746250 

-. 00 348750 

1.31 

1.32 

-.01524907 

-.01598805 

.28863737 
. 2966254 9 

.04042246 

.04351317 

-. 00381077 
-.00415061 

1.33 

1.34 

-.01672840 

-.01746909 

. 30449970 
. 31225793 

. 04673580 
.06009141 

-. 00 450710 
00 468025 


14 



TABLE A-2.— Coefficients for four-point integration — Continued 


1.35 

1.36 

“o0 1820911 
- .0 1894752 

. 31969018 
o 3 274185 6 

• 05 35 8099 
. 05720544 

“o 0052700 5 
-.00567648 

1.37 
1. 38 

- .0 1968337 
- . 0 20 41 575 

.33481727 

.34209259 

.06096556 
. 06486207 

".00609947 
" « Ob 653891 

1.39 

".02114380 

• 34924290 

. 06889560 

-.00699470 

1.40 

- . 02186667 

• 35626667 

. 07306667 

-.00746667 

1.41 

1.42 

".0 2258353 
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TABLE A-2. — Coefficients for four- point integration — Continued 
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Table A-3. Coefficients for six-point interpolation — Continued 
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TABLE A -3. — Coefficients for six-point interpolation — Continued 
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TABLE A-4. — Coefficients for six-point differentiation 
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TABLE A-4. — Coefficients for six-point differentiation — Continued 
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TABLE A-4. — Coefficients for six-point differentiation — Continued 
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TABLE A -4. — Coefficients for six-pomt differentiation — Continued 


<JO 

ro 

cn 

O 

CP 

M) 

ro 

CP 

in 

4 

<M 

O' 

o 

o 

CD 

K 

CM 

o 

CO 

M3 

CP 

O 

rH 

H 

CM 

rH 

lf\ 

on 

tH 

4 

CM 

CM 

CM 

ro 

ro 

CD 

O 

a 

CD 

CD 

« 

0 

c 

» 

A 

1 

1 

1 

• 

1 


in 

N* 

CO 

K 

in 

4 

M3 

r- 

rH 

O' 

N 

CM 

CO 

rH 

ro 

M3 

4 

iH 

M3 

ro 

CP 

CD 

ro 


ro 

4 

M3 

n- 

00 

o 

P- 

<P 

tH 

ro 

M3 

■H 

• 

T~t 

• 

CM 

a 

CM 

a 

CM 

• 


o 

CD 

o p- 

4 

rH 

CD 


4 

CD 

M3 

00 

ro 

CM 

4 

in 

rH 

CP 

in 

^ r 

CM 

in 

CM 

M> 

4 

eC 

rH 

4 

M3 

00 

4 

W N. 

ro 

00 

CO 

CO 


N- 


ro 

uD 

ro 

CP 

o 

ro 

vH 

CO 

CP 

CD 

ro 

rH 

in 

CQ 

CM 

ro 

CP 

CM 

CM 

(P 

ro 

ro 

4 

4" 

ro 

p- 

CD 

ro 

M) 

CP 

ro 

4 

4 

4 

4 

o 

o 

o 

CD 

o 

a 

a 

e 

Q 

a 

1 

P 

0 

♦ 

1 


O 

in 

h- 

0D 

N- 

o 

CP 

rH 

P- 

M3 

o 


N. 

P- 

LD 

o 

CM 

in 

M3 

rl 

o 

r-- 

4' 

rH 

00 

CM 

ro 

in 

N- 

00 

0G 

CD 

CM 

4 

M) 

CM 

ro 

ro 

ro 

ro 

a 

0 

0 

• 

9 


N. 

CD 

O 

o 

p- 

M3 

rH 

O 

rH 

M) 

Mj 

rH 

ro 

CM 

o 

M) 

rH 

M) 

CD 

rH 

M) 

4 

CP 

4 

CO 

C0 

CD 

■H 

ro 

4 

CP 

M3 

CM 

oc 

4 

M3 

M3 

MO 

in 

in 


CP 

ro 

M) 

ro 

CP 

<T 

00 

rH 

ro 

4 

in 

rH 

CD 

rH 

in 

o 

v D 

in 

M) 

00 

ro 

rH 

<T> 

M) 

CM 

CM 

in 

p- 

o 

ro 

in 

in 

in 

M3 

M3 

o 

o 

o 

o 

o 

0 

0 

a 

e 

0 

1 

P 

e 

1 

1 


in 

o 

in 

N. 

CO 

4 

in 

CP 

M3 

CM 

CM 

00 

CM 

co 

P~ 

4 

rH 

CP 

rH 

4 

ro 

P- 

00 

CO 

4 

o 

rH 

CM 

ro 

4 

cr 

tH 

ro 

in 

P- 

ro 

4 

4 

4 

4 

0 

• 

« 

0 

• 


4 

N. 

O 

o 

CD 

4 

M3 

rH 

o 

rH 

ro 

N 

ro 

CM 

CP 


N- 

rH 

N 

4 

(M 

CO 

P- 

CO 

4 

M) 

N 

CP 

rH 

4 

CD 

M3 

CM 

CP 

in 

in 

4 

4 

ro 

ro 


o 

CP 

ro 

M3 

ro 

in 

on 

ro 

M3 

co 

CM 

vH 

ro 

in 

N- 

rH 

ro 

ro 

o 

ro 

CO 

CM 

in 

K 

P~ 

in 

CO 

o 

CM 

4 

M3 

M3 

fv 

P- 

N. 

CD 

O 

CD 

O 

o 

a 

Q 

e 

« 

0 

0 

i 

0 

I 

» 

p- 

in 

-D 

in 

r- 

rH 

4 

CD 

4 

rH 

CP 

ro 

CO 

O' 

ro 

CM 


4 

N 

LA 

N. 

vO 

o 

CP 

ro 

4 

4 

4 

CM 

rH 

CP 

rH 

ro 

in 

N- 

4 

in 

in 

in 

in 

0 

a 

e 

a 

0 


N- 

4 

P- 

O 

O 

M3 

4 

vC 

rH 

CD 

rH 

CP 

4 

CM 

CP 

4 

4 

N 

CM 

CP 

in 

CM 

M3 

cr 

o 

N. 

rH 

in 

CD 

p- 

rH 

CO 

4 

rH 


ro 

CM 

CM 

CM 

tH 


CP 

o 

<P 

ro 

M3 

4 

o 

4 

GO 

M) 

00 

M3 

oo 

ro 

CP 

rH 

ro 

N 

ro 

00 

M3 

ro 

eO 

CM 

to 

M3 

oo 

CP 

rt 

CM 

P- 

p- 

P- 

00 

00 

CD 

CD 

CD 

CD 

CD 

0 


0 

o 

« 

1 

0 

0 

0 

I 


CO 

N 

in 

o 

in 

CM 

M3 

CP 

m 

4 

ro 

CM 

CM 

4 

M3 

rH 

CD 

M> 

ro 

in 

rH 

CM 

in 

rH 

oo 

CP 

M3 

CM 

CO 

CM 

CD 

CD 

CM 

ro 

tn 

in 

<.0 

M) 

M> 

M3 

a 

0 

0 

0 

a 


CD 


4 

P-. 

o 

rl 

M3 

4 

M3 

rl 

in 

CM 

M3 

ro 

4 

rH 

CO 

rH 

ro 

in 

ro 

M3 

ro 

ro 

OO 

4 

CM 

CM 

ro 

in 

4 

rH 

00 

in 

CM 

rH 

rH 

O 

o 

CD 


4 

rv 

Q 

CD 

o 

P- 

4 

P- 

O 

CD 

O 

N- 

4 

N. 

o 

O 

CD 

r- 

4 

P- 

O 

CD 

CD 

N- 

4 

4 

M3 

rH 

a 

rl 

M3 

4 

M3 

rH 

CD 

rH 

M3 

4 

M3 

rH 

o 

rH 

M3 

4 

M3 

rl 

CD 

rH 

MD 1 

4 

00 

CD 

P- 

ro 

M3 

<L> P*- 

ro 

ro 

M) 

in 

in 

in 

4 

M3 

in cp 

00 

M3 

P- 

O 

4 

CM 

CP 

4 

CP 

ro 

in 

4 

in 

M3 4 

N 

CM 

P- 

rH 

CM 

(p 

rH 

P- 

p- 

CD 

P^ 

00 

P'J 

4 

o 

4 

M3 

CD 

ro 

N. 

4 

M) 

ro 

M3 

M3 

ro 

CP 

ro 

GO 

ro 

CP 

CP 

rH 

cO 

rH 

CP 

in 

CD 

4 

CP 

in 

4 

OO 

ro 

O 

p- 

ro 

CP 

4 CP 

4 

03 

IO 

P~ 

CM 

M3 

rH 

P- 

CM 

CP 

in 

ro 

CM 

H 

rH 

ro 

M3 

CD 

ro 

ri 

00 

M) 

ro 

rH 

GO 

M3 

ro 

■H 

CO 

M3 

ro 

rl 

CC 

M3 ro 

rH (P 

P- 

in 

IO 

rH 

cn 

00 

CM 

CM 

vH 

rH 

rH 

rH 

O 

CD 

CD 

CD 

CP 

CP 

CP 

CP 

c0 

GO | 

CO 

CO 

p- 

P- 

n- 

P- 

p- 

MD 

M3 

a 

0 

a 

0 

0 

• 

0 

0 

• 

0 

» 

0 

0 

0 

o 

a 

« 

9 

« 

• 

o 

o 

« 

• 

0 

rH 

rH 

rH 

rl 

rH 

rH 

rl 

rH 

rH 

H 

1 

1 

1 

1 

1 

1 

t 

1 

I 

1 

i 

1 

i 

t 

1 

» 

1 

1 

1 

I 

1 

I 

1 

1 

1 
















tn 

ro 

CM 

ro 

in 

O 

in 

ro 

CM 

ro 

in 

o 

in 

ro 

CM 

ro 

in 

CD 

in 

ro 

CM 

ro 

in 

o 

m 

CD 

ro 

P- 

00 

in 

CD 

in 

co 

N- 

ro 

o 

in 

in 

IO 

CM 

00 

CD 

O 

o 

oO 

CM 

ro 

in 

in 

CD 

CD 

CP 

CP 

CM 

rH 

O 


00 

ro 

s- 

o 

ro 

o 

in 

fO 

o 

a) 

CO 

p- 

4 

4 

CP 

in 

CP 'CP 


M3 

o0 

M3 

CM 

o 

“ro 

m. 

rl 

ro 

P- 

in 

ro 

4 

4 

N 

N* 

CD 

rH 

in 

r- 

CVJ 

p- 

M3 

M3 

4 

M3 

in 

CM 

P~ 

i D 

rH 

rH 

rH 

CD 

Q' 

CP 

CD 

CM 

M) 

<M 

rl 

4 

CD 

CD 

in 

M3 

CM 

in 

MV 

ro 

P- 

rH 

in 

OC 

CM IT 

cO 

rl 

4 

M> 

CP 

ro 

M3 

CP 

ro 

N* 

rH M) 

rl 

M3 

CM 

CP 

M3 

4 

CM 

rH 

rH 

o 

CP 

CP 

00 

P- 

P- 

M3 

U' 

4 

4 

ro 

CM 

CM 

rl 

rl i 

CD 

CD 

<P 

CP 

OC 

OO ,00 

CM 

• 

CM 

0 

CM 

0 

CM 

« 

rH 

a 

rl 

0 

rH 

: 0 

i 

H 

V 

rH 

0 

rH 

0 

rH 

0 

rl 

0 

rH 

0 

rH 

0 

rl 

0 

rH 

0 

rH 

• 

rH 

• 

rH 

0 

rH 

« 

CD 

0 

CD 

• 

CD 

0 

o , 

0 > 

CD 

0 

rH 

O 

rH 

p- 

4 

P» 

t 

,tH 

O 

rH N, 

'4 

Ps. 

rl Olrl 

> N 

L J 

t 

CD 

rl 

P* 

4 

r-* 

rH 

in o 

fin 

tH 

ro 

iO 

in 

in 

M3; 00 

rH 

O 

O ICD 

rH 

or 

M3 


in 

CD 

M3 

ro 

rH fin 

M> cO UO 

4ir4 

vO'S- 

CPlrH 

o 

4 

CM 

rl 

CD 

'GO i 

i 4 j P- 

0° 

M3 

rH 

4 

4 

4 

4 tn 

CD 

o M3 

CP'CP 

M3 

CM 

^1 

fO CPi 

'M3 

<0 

CP M> N. 

1 in: CP 

rH ! 

CM 

ro 

4 

P- 

ro 

ro <o 

in CM cO 

4 fO 

M3 i 

CM 

N-i 

ro 


4 

LD 

M3 CO O I 

I 00 

M> 

M3| 

I M) 

N. 

CP 

CM 

p- 

ro io 

P- 

M3 

,4 


,CM 

CD 

CP 

p- 

M3 

4 


CM 

o CP 

co 

; MPin 

4 

1 no 

CVJ 

rl 

rH 

CD 

CD 

CD 

CM 

CM (CM CM 'CM 

(M 

rH 

rH, 

*-» 

rH 

h 

rllrlQ 

CD 

O 

o 

CD 

1 Q 

CD 

O 

O 

CD 

a 

O 

o 

0 

<=>1 

0 1 

o 

0 

O I 

CD 1 

CD 

CD 

CD 

iO 


!o 

t A 

CD .O 

a 

O j 

Oi 

CD 

O 

o 

CD 

CD 

CD 

O 

o 

O 

1 

1 

1 

'*! 

1 .* 

i 

1 

1 

1 

1 

1 

1' 

1 j 

1 

1 

1 

e 

1 

l 

1 

'' 

1 

1 

I 

• 

1 

♦ 

1 

0 

1 

1 

in 

M> 


ooj 

1 

CP 

O 

rH 

CM 

ro 

4 

lm 

M3| 

p- 

00 

CP 

1 

o 

rH 

CJ 

ro 

4 

in 

M) 

P^ 

oo 

Icr 

N 

r- 

N 


N~ 

00 

CO 

CO 

00 

00 

100 

GOj 

00 

CO 

oo 

<p 

cr 

CP 

a> 

CP 

cn 

CP 

CP 

CP 

kp 

0 

c 

I * 

• 

• 

0 

• 

0 

• 

• 

1 * 

• 

0 

0 

0 

o 

0 

0 

• 

9 

0 

• 

• 

0 1 

r. 

CM 

CM 

|CM 

CM 

(M 

CVJ 

CM 

CM 

ICM 

CM 

ICM 

CM i 

(M 

CM 

CM 

CM 

CVJ 

CVJ 

CM 

CM 

CM 

CM 

CVJ 

cm! 

irvj 


24 


OGQGQQOG .08333333 -.66666667 O.OCOQOGOQ .66666667 -.08333333 



EXAMPLES USING TABLES A-l AND A-2 


EXAMPLES USING TABLES A-3 AND A-4 


xn tuple 1.- Interpolate a value of Y at x=0.07 

von t lie table of base points 


r 

11 Yi 

-0.05 

1.00 

.00 

1.15 

.05 

1.25 

.10 

1 .25 

15 

1.10 

.20 

.75 

.25 

.15 

■=l + (0 07- 

-0 05)^0.10-0 05)=1 -1 


om table A-l: 


(C0= 

-0.072)xl.l5= 

-0.0828 

(Cl = 

.696)xl.25= 

.8700 

(C2= 

.424)xl.25= 

.5300 

(C3= 

~.048)xl.l0= 

-.0518 


y o 07 

1.2644 


Example 3. — Given the same base points as in 
example 1 interpolate at x=0.07. 

From table A-3: 


(C0= 

0.0144)x 1.00= 

0.0144 

(Cl — 

-.11361x1.15= 

-.1306 

(C2= 

.72641x1.25= 

.9080 

IC3= 

43841x1.25= 

.5480 

(C4 = 

-.07361x1.10 = 

-.0810 

( C5= 

.0080) x .75= 

.0060 


1 II 117 “ 

1 .2648 


Example *1. — Interpolate for the differential of y 
at .(=0.07 

From table A-4: 


U"0= 

0 01 800 lx 

1 00= 

0.0180 

tCl = 

,09733)x 

1.15= 

.1119 

l ("1— 

l.32267)x 

1.25 = 

1 .6533 

tr.i= 

1. 46400) x 

1 .25= 

1 .8300 

(C4 = 

.25933)x 

1.10= 

-.2853 

(C5= 

. 03867) x 

.75= 

.0290 




.0143 


ample 2. — Compute the integral from 2=0.00 
.r=0.07 


itltiMt > (HJ7 =(001 13 -’0 051=0.280 


<>m 0.00 to 0.05 
om table A-2: 


NOTES FC 

(f'0= -0.04167)xl.00= 

-0.0417 


(Cl= .54167)xl.l5= 

.6229 

TAH(l) 

(C’2= .541671x1.25= 

.6771 

(C3= -.041671x1.25= 

-.0521 



1.2062 

YINTdl 

om 0 05 to 0.07 


1 A1NT(1) 
YIMNT(I) 

om table A-2: 


INTFOHt 

(C0= -0.021871x1.15= 

-0.0252 

INTSIXt 

(Cl= .356271x1.25= 

.4453 


(C2= .073071x1.25= 

.0913 

nil 

(C3= -.007471x1.10 = 

-.0082 

K’OFFlI) 


.5032 

YAO(I) 

C 0 07 

J 0 Q0 =11.20(32+0 50321x0.05=0. 085-17 

X III) 


NOTES FOR PROGRAMS A-l, A-2, A-3, AND A-4 


Input of tabular values, or base 
points, on which interpolation is 
based . 

The interpolated values. 

Interpolated integral. 

Interpolated derivative. 

) Subroutine for four-point 
interpolation. 

1 Subroutine for six-point 
interpolation. 

Coefficients of equations 37 and 43. 

The matrices of equations 37 and 
43. 

Coefficients of equation 42. 

The index z of tables A-l, A-2, A-3, 
and A-4. 
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PROGRAM A-l . — Four-point interpolation and integration 


DIMENSION YINT (200 ) , Y AINT (200) 

COMMON/XX/TAP (100 ) 

READ<5»1000) (TAt) (I) ,1 = 3,22) 

1 0 0 0 FORMAT (20F4.2) 

T A P ( 1 ) = 0 • 1 
TAP (2) = 0. 1 
X = 0 .3 
MS = ~i 

no iooi T=i,i5i 
X=X+0.J1 

CALL INTFOR (X,Y ,YA , MS ) 

YINTCH =Y 
YAINT(I) = YA 
1001 CONTINUE 

WRITE <6, 100?) (I, YI NT (I), 1=1, 151) 

1 002 c ORMAT ( * INTERPOLATED FUNCT ION* / (8 ( 15, F 10 . 4> > ) 
WRITE (6, 1003) (I, YAINTCI) ,T = 1, 151) 

1 0 0 3 FORMAT <* INTERPOLATED INTEGRAL* / (fl < 15, F 10 . 4 )> ) 
STOP 
END 


SUBROUTINE INTFOR ( X, Y , Y A , MS) 

DIMENSION C(4) , ICOEF (4,4) 

C OMMON/ XX / T AP < 100 ) 

OATA ICOEF/4, -A, 5, -1,-6, 19, “14, 3, 6 , -16, 13, -3, -2, 5, -4,1/ 
AX=((X-TAB(1))/TAR<2) > 

M=I NT (AX) 

FX=AX-FLOAT(M>+l.Q 
IF(M.EQ.MS) GO TO 102 

no loo 1=1,4 
c m = o.o 
no ioo j= i , 4 

10 0 c (I) =C(I) + ICOEF (I , J) *TAB (M+JH) 

CA = C(i>+C (2) /2*0 + C<3)/3. 04-0(45/4.0 
MS=M 

102 Y=C<4) 

Y A= C ( 4 ) /4. 0 
HO 101 K=1 ,3 

Y = Y * F X 4- C ( 4 “ K ) 

101 YA = YA*'=X4-C (4-K) /(4-K) 

Y=Y/2. a 

Y A= (YA*FX-CA)*T AP( 2) /2.0 
RETURN 

END 
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PROGRAM A-2. — Generation of four-point coefficients 


1 1MENSt ON YINf <~10'l ,4)~, YAINT (101,4) ,IOOEF(4,4) , XI (101) , YA0(4) 
. . 0 A TA .Ip _0E F _/4 j i ” 8j 5, -1,-6, 19,-14,3, 6, ~16, 13, “3, -2, 5, -4,1/ 

YAO(l) =-17.0/24.3 
Y A0(2)_=5. 0/24.0 
YAO (3) =-19.0/24.0 " 

YAO (4) =7. 0/24.0 

00 10GD 'j = l,l6f 

X = FLOA T { J ) /100. JfC .99 

XI (J) =X 

OO 1001 1=1,4 

Y INT ( J , I) = ICOEF (4 , I) ‘ 

OO 1C u 2 K = l,3 

10J 2 Y'INT <j,I) =YINT( J,I j*X+IC0EF (4-K, I) 

YINT(J,I) =YINT<_J,I )/2.0 
YAINT (J,I)=ICOEF (4, I ) /4 . 0 
DO 10 0 5 K = 1 , 3 

10G3 Y AINT ( y , I ) = ( YAINT ( J, I ) *X) + ( FLOAT (ICOEF ( 4-K ,1 ) ) / (4-K) ) 

YAINT (J, I )_= (YAINT ( J,I)/2.0> *X+YAO(I) 

10 01 ‘CONTINUE 
1 00 0 O ONTIN'JE 

OO 1003 J=l,101 

IF( J.EQ.l .OR. J.EQ. 51) GO TO 1006 
3 O TO 1007 
1006 WRITE (6, 1C 08) 

1008 FORMAT (*1 ¥ ) 

1 0 u 7 IF (MOO (J, 5_). EQ. 1J GO TO 1009 
GO TO 101C 

1009 NRITE(5, iCill) 

1011 FORMAT <*0*) 

1010 WRITE(6,JU12) ( X_r < J > , (YINT(J,I) ,1 = 1,4)) 

1012 FORMATC* '* \ ( 35 X , F£> i 2~, 4 FI 2 . 8 ) ) 

1006 C ONTINUE_ 

DO 1020 J =1,101 

IF(Jj,E3.1 .J3R.J.EQ.51) GO TO 1013 
GO TO" 101 4 

1013 WRITEJ 6, 1 0C8) 

1014 1 FfMOCft 3V 5> • tn.“£> GO TO 1016 
GO TO 1017_ _ 

1016 WRIT? (6,1 011) 

1017 WRITE (6,1012) (X_I<J) , ( YAINT ( J, I ) , 1 = 1 ,4) ) 

10 20 CONTINUE" 

STOP 

“ 'END ‘ 
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PROGRAM A-3. — Six-point interpolation and differentiation 


DIMENSION YINT (200 > YPINT (200) 

COMM ON /XX/TAP1 1-00) 

REA0<5, 10 00> <TA 9(1 > ,1 = 3,22) 

1000 FORMAT { 20 F4. 2) 

T A D ( 1 ) = 0 , 1 
TAP<2)=0.1 

X = 0 . 3 
M S= - 1 

00 tool 1=1,151 

X=XfO .31 

CALL INTSIX(X»Y ,YP ,MS) 

Y INT ( I ) =Y 
YPINT(I)=YP 

1001 CONTINUE 

WRITE (6, 1002X1, YI NT (I), 1=1, 1*1) 

10 0 2 F ORMAT <* INTERPOLATED FUNCTION */< B ( 15, F 10 . 4) ) ) 
WRITE (6, 100 3) (I, Y PINT (I) ,1=1,151) 

1003 FORMAT (* INTERPOLATED DE RI V AT I VE V < 8 <15 , FiQ . 4 )1 ) 
STOP 
END 


SUBROUTINE INTSIX(X,Y , YP , MS) 

DIMENSION C ( 6) ,IC0EF (6,6) 

CCMM0N/XX/TAP(130 > 

DATA ICOEF/432, -918,765, -313, 63, -5, -20 4 0,4436,-3754, 15 5 1,-31 4, 25, 

1- * 080, -875 2, 7414, -3 07ft, 62 6, -50 ,-4030,0712,-7356, 305 8,-624,50,2 040, 

2 - 4346 , 3661 ,- 1521 , 311 ,- 25 ,- 408 , 868 ,- 730 , 303 ,- 62 , 5 / 

AX=((X-TAB<i) ) /TAP (2) ) 

M=INT( AX) 

F X= AX-FLOAT (M) +2. 0 
IFCM.EQ.MS) CO TO 102 
DO 100 1=1,6 

c m=c.o 

DO 100 J= 1 , 6 

100 C C I > =C f I) +IC0EF (I , J) *TA9(M+J) 

MS = M 

10? Y=C (6) *F<+C(5> 

YP=5 • 0 *C ( 6) 

DO 101 K= 1 , 4 
Y=Y*FX+C(5-K) 

101 YP=Y d *FX+ (5-K)*C(6-K) 

Y=Y/24, 0 

YP = YP/2 4. 0 /TAB (2 ) 

RETURN 

END 
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PROGRAM A-4. — Generation of six-point coefficients 


DIMENSION YINT (101*6) ,YPINT(ioi,6) ,ICOF.F<6,6),XI(101) 

DATA ICQEF/4 32, -91 8,7 65, -313, 63, -5 , -20 40, 44 36, -375 4, 1551, -314,25, 
14 0 8 C,-? 7 5 2, 7414 , -3 07 8,626 , -50 , -40 BP, 87 12, -7356, 305 8, -624, 50, 20 40 , 
2- 43 46, 3 66 1 ,-152 1,3 11, -25, -400,86 8, “730, 303, -62, 5/ 

00 100 1 J = 1 ,101 
X = F.LOAT (J ) 7-lflQa 0tl,99 
X I ( J ) =X 

DO 1001 1 = 1,6 . .. 

Y INT (J , I) =TCOEF (6 , I) 

D Q- 10 0 2 K = 1 , 5 . . . 

L002 YINT (J, I) =YINT< J , I > *X+ICOEF (6-K , I) 

YINT (J,I) =YINT( J,I)/24.0 
YPINT ( J,I ) = ICOE F (6 ,1) *5. 0 
10 1003 K = 1 , 4 

1 0 0 3 YPINT (J ,1 ) =Y n I NT ( J , I ) *X+TCOEF ( 6-K, I ) * ( 5-K ) 

YPINT U,I) =Y^INT (J ,1) /24.0 
1001 DONTIIMIJ^ 

L 0 Q 0 CONTINUE 

DO 1005 J=l,101 

IF(J.ET,1. OR. J. EQ. 51 ) GO TO 1006 
CO Tp 100 7 
1006 WRITE(5,100a> 

1008 FORMAT ( *1 * ) 

L 0 0 7 tF(MOD(J,5)»EQ»l) GO TO 1009 
CO TO 10 1C 
tOOD WRITE ( 6 , 1 0 11) 

LOii F CRMA T ( *0 * ) 

1010 WRITE (6, 10 12) (XI (J) , ( YINT ( J , I ) ,1=1,6)) 

L 0 1 2 FORMATS " , (35X ,F5 .2, 6F12. 0) ) 

L 0 0 n CONTINUE 

DO 1020 J= 1, 1 0 1 

I F ( J.F1.1 .OR.J.EO.51) GO TO 1013 
CO TO 1014 
I 013 WRITE (6,10 08) 

L 0 1 4 IF(MO0(J,5) .FO. 1) CO TO 1016 
CO TO 1017 
1016 WRITE (6,1011) 

L 0 1 7 WRITE (6 ,1012) (XI (J) , (YPINT (J,I ) ,1 = 1 ,6) ) 

1020 CONTINUE 

STOP 
c Nfi 
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APPENDIX B. — DATA ARRANGEMENT AND COMPUTER PROGRAMS FOR 

SMOOTHING 


TAHI.E B-l. — Data arrangement for two-point, smoothing 


Measured Kase-point Scale Values of coefficients 1 

<lata points positions vanate for base points 


X 

V 

111 X 

z 

YB 1 

YB2 

YB3 

YB4 

V 



6.4 







0 412 

1.50 


0.0G 

0.94 

0.00 



1.50 

4 Hi 

1.09 


.08 

.92 

.08 



1.69 

458 

2.02 


.29 

71 

.29 



2.02 

,450 

2 20 


25 

.75 

.25 



2.20 

.404 

2.50 


.32 

.08 

.32 



2.50 

.500 

2.50 


.50 

• .50 

.50 



2.50 

500 

2.70 


.50 

.50 

.50 



2.70 

.514 

2.89 


.57 

.42 

.57 



2.89 

.560 

3 25 


.S3 

.17 

.83 



3.25 

.574 

3 58 


,87 

.13 

,87 



3.58 

.000 

3.02 

6 

1.00 


1 00 



3.62 

.024 

3 60 


.12 


.88 

0.12 


3.66 

.020 

3.85 


13 


.87 

.13 


3.85 

.000 

3.55 


.30 


.70 

.30 


3,55 

000 

4,04 


.30 


70 

.30 


4.04 

,080 

3 88 


,43 


.57 

43 


3.88 

.098 

3.60 


.49 


.51 

.49 


3.60 

.710 

3.98 


58 


.42 

.58 


3.98 

750 

3.70 


75 


.25 

.75 


3.70 

774 

3.87 

,8 

.87 


.13 

.87 


3.87 

.800 

3 58 


.03 



.97 

,03 

3.58 

.844 

3.09 


.22 



.78 

22 

3.69 

.800 

3 29 


30 



.70 

.30 

3.29 

.800 

3.50 


.33 



.67 

.33 

3 56 

.880 

3 27 


.43 



.57 

.43 

3.27 

.910 

3.10 


.55 



.45 

.55 

3.10 

,914 

3 23 


.57 



.43 

.57 

3.23 

.950 

2,74 


.78 



.22 

78 

2,74 

.900 

2.98 


,80 



.20 

,80 

2,98 

.980 

2 55 

1.0 

.93 



.07 

.93 

2.55 





X\ 

*2 

*3 

*4 

y 


Conventional regression symbols 


1 liquations 12 and 13. 
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TABLE B-2. —Data airangement for four-point smoothing 


Base* point Scale Values of coefficients 1 

positions variate for base points 

in x £ YB1 YB2 TO3 YB\ YB5 YBb Y 

0 A 

.6 


1.06 

-0.0265 

0.9913 

0.0369 

-0.0017 



1.50 

1.08 

-.0339 

.9848 

.0520 

-.0029 



1.69 

1.29 

-.0731 

,8263 

.2766 

-.0299 



2.02 

1.25 

-.0703 

8672 

.2266 

-.0234 



2 20 

1.32 

-.0740 

.7932 

.3156 

-.0348 



2.50 

1.50 

-.0625 

5625 

.5625 

-.0625 



2.50 

1.50 

-.0625 

5625 

.5625 

-.0625 



2.70 

1.57 

-.0527 

.4655 

.6570 

-.0699 



2.89 

1.83 

-.0120 

.1354 

,9351 

-.0586 



3.25 

1 87 

- 0074 

.0955 

.9610 

- 0492 



3.58 

2.00 

.0000 

.0000 

1.0000 

.0000 

0.0000 


3.62 

1.12 


-.0465 

.9666 

.0862 

-.0063 


3.66 

1.13 


-.0492 

.9610 

.0955 

-.0074 


3.85 

1.30 


-.0736 

.8155 

.2895 

-.0315 


3.55 

1.30 


- 0735 

,8155 

.2895 

- 0315 


4.04 

1.43 


-.0699 

.6570 

.4655 

-.0527 


3.88 

1.49 


-.0637 

.5762 

.5487 

-.0612 


3.60 

1.58 


-.0512 

.4517 

.6701 

-,0706 


3.98 

1.75 


-.0234 

.2266 

.8672 

-.0703 


3.70 

1.87 


-.0074 

.0955 

,9610 

-.0492 


3,87 

1.03 


.0000 

-.0141 

.9978 

0168 

-.0004 

3.58 

1.22 



-.0669 

.8950 

.1908 

-.0189 

3.69 

1.30 



-.0735 

.8155 

.2895 

-.0315 

3.29 

1.33 



-.0741 

.7417 

.3289 

-.0365 

3,56 

1.43 



-.0699 

.6570 

.4655 

-.0527 

3.27 

1.55 



-.0557 

.4933 

.6304 

-.0681 

3.10 

1.57 



-.0527 

.4655 

.6570 

-.0699 

3,23 

1.78 



-.0189 

.1908 

.8950 

-.0669 

2.74 

1.80 



-.0160 

.1680 

.9126 

-.0640 

2.98 

1.93 



-.0023 

.0443 

.9883 

-.0303 

2.55 



*2 

*3 

X 4 

*5 


Y 


Conventional regression symbols 


i 

p>om table A l . 
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Table B-3. —Data arrangement for six-point smoothing 


— 

— 




Values of coefficients 1 




Base point 

Scale 





foi base points 




_____ 

positions 
in r 

variate 

YB 1 

ym 

YB3 

YBA 

YB3 

YBS 

YB1 

Yim 

y 

04 

6 











8 

2.06 

0.0048 

-0 0373 

0.9949 

0.0429 

-0,0054 

0.0001 



1 50 


2 08 

0062 

- 0483 

9907 

.0588 

-.0075 

.0001 



1 69 


2 29 

.01-19 

-.1144 

.8566 

.2842 

-.0454 

.0040 



2 02 


2 25 

,0143 

- 1090 

.8948 

2332 

- 0360 

0028 



2 20 


2,32 

0151 

-.1163 

,8245 

.3245 

-.0528 

0050 



2 50 


2 50 

.0117 

- 0977 

.5859 

.5859 

-.0977 

.0117 



2 50 


2 50 

.0117 

-.0977 

5859 

.5859 

- 0977 

.0117 



2 .70 


2 57 

.0092 

- 0813 

4826 

6858 

-.1101 

.0138 



2.89 


2.83 

0010 

-.0202 

1421 

.9536 

-.0882 

0115 



3 23 


2 87 

.0005 

-.0139 

1026 

.9739 

- 0726 

.0095 



3 58 

1.0 

3 00 

0000 

.0000 

.0000 

1.0000 

0000 

ooou 



3.62 

2 12 


.0089 

-.0681 

9780 

0933 

- 0125 

0 0004 


3.66 


2.13 


0095 

- 0726 

.9739 

1026 

- 0139 

0005 


3 85 


2.30 


0150 

-.1152 

8462 

.2975 

-.0478 

.0043 


3 55 


2,30 


.0150 

-.1152 

.8462 

.2975 

- 0478 

.0043 


4 04 


2 43 


0138 

-.1101 

.6858 

.4826 

-.0813 

0092 


3.88 


2 49 


.0121 

- 0997 

.6005 

.5713 

-.0955 

.0114 


3.60 


2.58 


0088 

-.0787 

.4678 

.6995 

-.1114 

.0140 


3.98 


2 75 


.0028 

- 0360 

.2332 

8948 

-.1090 

.0143 


3.70 


2.87 


0005 

- 0139 

.1026 

.9739 

- 0726 

,0095 


3 87 

1 2 

2 03 



0025 

- 0194 

.9988 

.0207 

-.0026 

0.0000 

3 58 


2.22 



.0135 

-.1028 

9197 

1972 

-.0296 

.0020 

3.69 


2.30 



0150 

-.1152 

.8462 

.2975 

-.0478 

.0043 

3.29 


2.33 



.0151 

- 1166 

.8132 

.3383 

- 0554 

,0054 

3.56 


2.43 



.0138 

-.1101 

.6858 

4826 

-.0813 

.0092 

3.27 


2 55 



.0099 

-.0862 

.5122 

.6579 

-.1071 

0133 

3,10 


2.57 



0092 

-.0813 

.4826 

.6858 

-.1101 

,0138 

3 23 


2.78 



.0020 

-.0296 

1972 

,9197 

-,1028 

0135 

2 74 


2.80 



.0016 

-.0256 

.1744 

9344 

-.0976 

.0128 

2 98 

1 4 

2.93 



0001 

-.0064 

.0508 

.9930 

-.0429 

.0055 

2,55 

1.6 

1.8 













*i 

*2 


*4 

*5 

*6 

*7 

*8 

Y 






Conventional regression 

symbols 





*From table A-3 


NOTES FOR PROGRAMS B-l AND B-2 

X(I) Input conventional re- 

gression data from 
tables B-2 and B-3. 

SSX(I) Sums of squares and 

products. 


AAU) 

MINV(AA,N,I>,L.M) 

HU) 

INTFOR(V,¥,YA,MS) 

INTSIX(V,Y,YP,MS) 


Matrix to be inverted. 
System library routine 
for matrix inversion. 
Regression coefficients 
or base points. 

See program A-l. 

See program A-3. 
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PROGRAM B-1. — Four-point smoothing 


DIMENSION X (30, 10) ,SSX (10 , 10) ,L(10) ,H(10) » B ( 1 0) ,YINT (60) ,YAINT ( 0 ) 
1 , Art ( 6 , G ) 

COMMON/XX/ TrtB(lOO) 

READ (5, 10 00) <(X(I,J) ,J=1,7> ,1=1,30) 

1000 FORMAT ( 7F5. 4) 

DO 1002 I =1,10 

no 1002 j=i,io 

loo? ssxd, J) = o.o 

DO 100 ( 1=1,7 
DO 1004 J = 1 , 7 
DO 1005 K = 1 , 30 

1005 SSX (I»J)=SSX(I,J)+X(K,1)*X(K,J) 

1004 CONTINUE 

1 u 0 3 CONTINUE 

DO 10 0G 1 = 1,7 
DO 100 6 J = 1 , 7 

1006 SSX (J,I) = SSX (I, J) 

WRITE (6,10 10 )( (SSX (I, J) ,1=1,7) ,J = 1, 7) 

1010 FORMATS *, (7F12.fi)) 

DO 1020 1=1,6 

DO 1020 J =1 , 6 
1020 AA(I,J)=SSX(I,J) 

CALL MINV (AA,6,Q,L ,M) 

WRITE(6,1011) ( < A A ( I , J ) ,1=1,6) »J = 1,6) 

1011 FORMAT (* *, (6F13.fi)> 

DO 1007 1=1,6 

D (I )=G . 0 
DO lOQfi J = 1 , 6 

1006 D (I ) =B ( I) +AA(I,J) ¥ SSX(J,7) 

1007 CONTINUE 

WRITE (6,1011) (B( I) ,1=1,6) 

TAR (1) =0. 4 
TAR (2) = 0. 2 
DO 1012 1=3, fi 

1012 T AB ( I ) = B ( 1-2) 

V = 0. 6 
MS=-1 

DO 101 i 1 = 1,31 

CALL INTFOR(V»Y*YA»MS) 

Y INT ( I ) =Y 
YAINT(I)=YA 
V=V+0. 12 

1013 CONTINUE 

Y AINT ( 1) = 0 » 0 

W RITE ( 6 ,1 0 14) (I, YI NT (I), 1 = 1, 31) 

1014 FORMAT (* INTERPOLATED FUNCTION*/ (8(15, F10. 6))) 

WRITE (6, 10 15) (I, YA INT (I), 1 = 1, 31) 

1015 F ORMAT ( * INTERPOLATED INTEGRAL*/(8(I5»F10*6))> 

STOP 

END 
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PROGRAM B-2. — Six-point smoothing 

DIMENSION X(30,lu),SSX(10,10),L<10) ,“M( i oTTsTl CITY! NT <60 ) ,YPINT (6 0 ) 
1,AA(8,3) 

COMMON/XX/ TAB(100 ) 

READ (5,10 0 0) ((X (I, J) ,J=1.,9> ,1 = 1,30) 

100 0 FORMAT (9F5. 4) 

90 1002 1=1,10 
' 90 1002 J= 1,10 
jL_0 0 2 SSX(I,J)=0.0 
90 1003 1=1,9 
DO 1004 J=I,9 
"90 1005 K=1 ,30 

1005 SSX(I,J)=SSX(I,J)+X(K,I)*X(K,J> 

100 4 CONTINUE “ ' 

100 3 CONTINUE 

90 1006 1=1,9 
90 1006 J=I ,9 

1006 SSX(J,I)=SSX(I,J) 

WRITE (6, 10 10) <<SSX(I, J) ,1=1,9) ,J=1,9) 

1010 F ORM AT ( * *,(9F12.8>> 

90 1020 1=1,8 

90 1020 J=i,fl 
1020 AA (I , J) =SSX (I, J) 

CALL MINV(AA,8,9,L,M) 

WRITF(6,10U) ( ( AA ( I , J> ,1 = 1,8) ,J = 1,8) 

1011 FORMAT** * , ( 8F1 3 , 8 ) ) ~ ' ‘ 

90 1007 1=1,8 

9 ( I > = 0 . 0 
DO 100 3 J = 1 , 8 

1 00 8 9 (I)=R(I) + A A ( I , J) * SSX ( J , 9 ) 

1007 CONTINUE 

WRITE(6, 1011) (B(I) ,1=1,8) ' ' ~ 

TAB(l) =0. 4 
TAR (2) = 0. 2 
00 1012 1=3,10 

1012 T AR(I) =B< 1-2) 

V=0 .8 

M S= -1 

90 1013 1=1,31 

CALL INTSIX(V,Y,YP*MS) 

YINT(I) =Y 
YPINT ( I ) = YP 
V=V+0.02 

1013 CONTINUE 

WRITE (6, 1014) (I,YINT (I) ,1=1,31) 

1014 FORMAT (* INTERPOLATED FUNCTIO N*/ < 8 ( 15 , F 10 . 6 > ) ) 

WRITF (6,1 015) (I, YPINT (I) ,1=1,31) 

10l£ FORMAT (* INTERPOLATED DERIVATIVE* /(8(I5»F10.6) ) ) 

STOP 

ENO 


* 1976 GPO 1750-S-671 -583/49 
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